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Abstract. We combine the idea of confluent drawings with Sugiyama 
style drawings, in order to reduce the edge crossings in the resultant 
drawings. Furthermore, it is easier to understand the structures of graphs 
from the mixed style drawings. The basic idea is to cover a layered graph 
by complete bipartite subgraphs (bicliques), then replace bicliques with 
tree-like structures. The biclique cover problem is reduced to a special 
edge coloring problem and solved by heuristic coloring algorithms. Our 
method can be extended to obtain multi-depth confluent layered draw- 
ings. 



1 Introduction 

Layered drawings visualize hierarchical graphs in a way such that vertices are ar- 
ranged in layers and edges are drawn as straight lines or curves connecting these 
layers. A common method was introduced by Sugiyama, Tagawa and Toda |25 | 
and by Carpano Several closely related methods were proposed later (see 

e.g. [niTiTnE2imiii-) 

Crossing reduction is one of the most important objectives in layered draw- 
ings. But it is well known that for two-layer graphs the straight-line crossing 
minimization problem is NP-complete ^1]. The problem remains NP-complete 
even when one layer is fixed. Jiinger and Mutzel |15) present exact algorithms 
for this problem, and perform experimental comparison of their results with 
various heuristic methods. Recently new methods related to crossing reduction 
(ESCEEaini) have been proposed. 

However when the given two-layer graph is dense, even in an optimum solu- 
tion, there are still a large number of crossings. Then the resulting straight-line 
drawing will be hard to read, since edge-crossing minimization is one of the most 
important aesthetic criteria for visualizing graphs [24 . This give us a motivation 
for exploring new approaches to reduce the crossings in a drawing other than 
the traditional methods. 

In addition, it is sometime of interest to find the bicliques between two layers. 
For example in the drawing of a call graph, it is interesting to find out which set 
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of modules are calling a common set of functions and what are those common 
functions. Call graphs are usually visualized as layered drawings. However it is 
hard to learn this information from layered drawings by traditional Sugiyania- 
style approaches, especially when the input graphs are dense. 

Our previous work 5 introduces the concept of confluent drawings. In |S] 
we talk about the confluent drawability of several classes of graphs and give a 
heuristic for finding confluent drawings of graphs with bounded arboricity. In 
this paper we experiment with an implementation of confluent drawings for the 
layered graphs. However we relax the constraint of planarity and allow crossings 
in the drawings, while it is not allowed to have crossings in a confluent drawing 
in our previous deflnitions. 

We are aware of the Edge Concentration method by Newbery |2I] . Edge Con- 
centration and our method share a same idea of covering by bicliques. But in 
Newbery's method, dummy nodes (edge concentrators) are explicit in the draw- 
ing and treated equally as original nodes, which causes the nodes' original levels 
to change. In our method dummy nodes are implicit in the curve representation 
of edges and the original levels are preserved. Furthermore, our method uses a 
very different algorithm to compute the biclique covers. 

2 Definitions 

In this section we give definitions for confluent layered drawings. The definitions 
almost remain the same as in our previous confluent drawing paper, except that 
the planarity constraints are dropped. Fig. Ogives an idea of confluent layered 
drawings. Edges in the drawing are represented as smooth curves. 




Fig. 1. An example confiuent layered drawing. 



A curve is locally monotone if it contains no sharp turns, that is, it contains 
no point with left and right tangents that forms a angle less than or equal to 90 
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degrees. Intuitively, a locally-monotone curve is like a single train track, which 
can make no sharp turns. Tracks are the union of locally-monotone curves. They 
are formed by merging edges together. 

A drawing A formed by a collection of tracks on the plane is called a confluent 
drawing for an undirected graph G if and only if 

— There is a one-to-one mapping between the vertices in G and A, so that, 
for each vertex v € V{G), there is a corresponding vertex v' £ A, and all 
vertices of G are assigned to distinct points in the plane. 

— There is an edge {vi,Vj) in E{G) if and only if there is a locally-monotone 
curve e' connecting v'^ and v'j in along tracks in A. 

The directed version of a confluent drawing is defined similarly, except that 
in such a drawing the locally-monotone curves are directed and in every track 
formed by the union of directed curves, the curves must be oriented consistently. 

Self loops and parallel edges of G are not allowed in our definitions, although 
multiple ways of realizing the same edge are allowed. Namely, for an edge in 
the original graph, there could be more than one locally monotone path in the 
drawing corresponding to this edge. 

We apply the idea of confluent drawings on layered graphs. Particularly, in 
the resultant confluent drawing, we replace bicliques in a biclique cover of a two- 
layer graph G = {U, L, E) by tree-like structures and draw them with smooth 
curves. As we can see in Fig. ^ our method can greatly reduce the crossings in 
the drawings of dense bipartite graphs. Additionally, nodes of a biclique can be 
easily identified by following the smooth curve paths. 

Since it is valid to have more than one confluent path between two nodes u 
and I in the confluent drawing when (u, I) G E, as defined above, it is straightfor- 
ward that a confluent layered drawing can be obtained by computing a biclique 
cover G of G, then visualizing each biclique in C as a tree-like structure. We 
show how to compute a biclique cover of G in the next section. 

3 Computing Biclique Cover of Bipartite Graphs 

Fishburn and Hammer 9 show that the biclique cover problem is equivalent to 
a restricted edge coloring problem. This coloring is not much useful for general 
graphs. However, it has a nice result for triangle-free graphs, and since bipartite 
graphs belong to the class of triangle-free graphs, an immediate result is that 
this type of edge coloring can be used to find a biclique cover of a bipartite 
graph. This result is useful in layered drawing because the edges between any 
two layers in such a drawing induce a bipartite subgraph. 

An edge coloring c: E <~ {1, 2, . . . , fc} for G = {V, E) is simply-restricted if 
no induced is monochromatic and the vertex-disjoint edges in an induced 
Pi or Cf have different colors. Fig. [2 shows the conditions that such induced 
subgraphs of a simply-restricted edge coloring must satisfy. 
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Fig. 2. The required conditions of induced subgraphs of a simply edge coloring. 

Let d{G) denote the bipartite dimension of G, which is the minimum cardi- 
nality of a biclique cover of G. Let Xs (G) be the chromatic number of a simply- 
restricted edge coloring of G. Xs{G) is if i? = 0; otherwise, it is the minimum k 
for which G has a simply-restricted coloring c: E ^ {1,2, . . . , k}. The following 
theorem states the equivalence of d{G) and Xs{G) for triangle-free graphs. 

Theorem 1 (Fishburn and Hammer f^) 

d{G) = Xs{G) for every triangle-free graph. 

Let Ej be the set of edges with color j in a simply-restricted edge coloring 
for a triangle-free graph G. As we can see in the second part of the proof of The- 
orem ^ (omitted here, included in the full version of this paper), Ej is included 
in the edge set of a biclique subgraph of G. Therefore, every edge set of a single 
color induces a biclique subgraph of G. By computing a simply restricted edge 
coloring we can get a biclique cover of G. 

Because it is known that the problem of Covering BY Complete Bipar- 
tite Subgraphs is NP-hard (Garey and Johnson ^ GT18), it is unlikely to 
have efficient optimization algorithms for finding the minimum biclique cover of 
a bipartite graph. Thus we only focus on fast heuristics for computing a near- 
optimal biclique cover. 

The simply-restricted edge coloring problem can be transformed into a vertex 
coloring problem. So, instead of devising a special algorithm for the simply- 
restricted edge coloring, we can choose to use one of the existing vertex coloring 
algorithms. Well known heuristic algorithms for vertex coloring include Recursive 
Largest First (RLF) algorithm of Leighton [TS|, DSATUR algorithm of Brelaz [2]. 
For more about heuristics on graph coloring, see Campers et al. [3]. 

The above method of computing a biclique cover by coloring doesn't distin- 
guish between two kinds of bicliques: Kp^i and Ki^r, where p,q,r > 1. So if 
we are more interested in finding out the set of common callers and callees, we 
would need to give higher priority to Kp q than Ki r when covering the edges. 

After the biclique cover of the two-layer bipartite graph is computed, each 
biclique in the cover is drawn as a tree-like structure in the final drawing. Doing 
this repeatedly for every two adjacent layers, we can get the drawings for multi- 
layer graphs. 

The time complexity of the algorithm depends on the coloring heuristic sub- 
routines. For a graph with a set of vertices V, both the RLF algorithm and the 
DSATUR algorithm run in worst case 0(|yp) time. There are some other faster 
coloring heuristics with 0(|yp) time, but their output qualities are worse. Sup- 
pose we have a two-layer bipartite graph G = {V, E). The transformation from 
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the simply-restricted edge coloring into vertex coloring version takes 0{\Ef) 
time. Using RLF or DSATUR costs 0{\E\^), thus the total time is 0{\E\^). 

4 Layout of the bicliques 

We described how to compute a biclique cover of a two-layer bipartite graph in 
the previous section. Now it is time to show how the bicliques are laid out. In 
the confluent layered drawings, each biclique in the biclique cover is visualized 
as a tree-like structure, as in Fig. ^ Now here are the questions. What are the 
best positions to place the centers of the tree-like structures? How to arrange 
the curves so that they form confluent tracks defined in Section [^P 

4.1 Barycenter method to place centers 

In the case where the positions of nodes in the upper level and lower level are 
fixed, one would like to put the center of a tree to the center of the nodes 
belonging to the corresponding biclique. For example, in Fig. |31 the drawing on 
the left is visually better than the drawing on the right. Firstly it has better 
angular resolution and better edge separation. Secondly it is easier for people to 
see the biclique as a whole. Then the next question is: what does the center of 
those nodes mean? In our method, the natural candidate position for a center 
of the tree-like structure is the barycenter, i.e., the average position, of all the 
nodes in this biclique. 




Fig. 3. Good-looking tree and bad-looking tree with centers placed differently. 

It looks bad too if these tree centers stay very close to each other. So we need 
to specify a minimum separation between two centers. 

The above requirements can be formulated into constraints: 
1. A tree center stays within the range of its leaves. 

min Xij < Xi < maxxy , 

j 3 

where Xi is the x-coordinate of the i**^ tree center Ci, and Xij is the x- 
coordinate of the j^^ leaf of a. 
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2. The distance between any two centers is greater than or equal to the mini- 
mum separation. 

V« ^ i, \xi - x.j\>5 
where 5 is some pre-specified minimum distance. 

Under these constraints, we want a tree center to stay as close as possible to 
the barycenter of all its leaf nodes, i.e., we want to minimize 'Yliii.^i 
subject to the above constraints. This is a Quadratic Programming problem, and 
unfortunately it is NP-hard (Garey and Johnson [151, MP2). 

Since it is unlikely to have efhcient algorithms for solving this optimization 
problem, and a small deviation of a tree center from the perfect position won't 
cause too much displeasure, we use instead a very simple heuristic method to 
place the tree centers. We first assign to each tree center the x-coordinate of the 
barycenter of its leaves. Then we sort tree centers by their x- value. The third step 
is to try to place these tree centers at their x-coordinates one by one. Assume 
we have k centers to place. Start from the j*'' center, where j = [^\. Place 
center j at its barycenter, then try to place centers one by one in the following 
order: j — 1, j — 2, . . . , 1. If constraint 2 is violated, the violating center is placed 
the minimum distance away from the previous placed center. Tree centers to the 
right of center j are placed similarly in the order of j + 1, j + 2, . . . , fc. It is easy to 
see that the running time of the barycenter method is dominated by the sorting 
of the tree centers. 

4.2 Placing tree centers to reduce crossings 

Alternatively, one might want to place these centers on positions such that the 
total number of edge crossings is as few as possible, especially in the case where 
nodes of upper level and lower level are not fixed. If this is the main concern, we 
can place the tree centers in another way in order to reduce the edge crossings. 

After the biclique cover of a two-layer graph G = {U,L,E) is computed, we 
construct a new three-layer graph G". We treat these tree centers as nodes of a 
middle layer. The set of vertices includes three levels: an upper layer U' = U, a 
middle layer M consisting of tree centers, and a lower layer V = L. The edges 
of G' are added as follows: for each biclique Bi in the biclique cover, add one 
edge between the tree center node and each node u <eU that belongs to Bi. 
Similarly add one edge between rm and each node I G L that belongs to Bi. 

Now a two-layer graph of the original problem is transformed into a three- 
layer graph G'. Straight-line crossing reduction algorithms can be applied on G' . 
After the crossing reduction, we obtain the ordering of nodes in each of the three 
layers. The orderings will be used to compute the positions of nodes and tree 
centers in the final drawings. Note that when crossing reduction method is used 
to place tree centers, it is not always true that a tree center always stays within 
the x-range of its leaves, i.e., bad centers hke the one in Fig.Elcould appear. 

Here we are using straight-line edge crossing reduction algorithms for our 
confluent layered drawings with curve edges. Readers may suspect the equality 
of the crossing number in the straight-line drawing for the new three-layer graph 
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G' and the crossing number of our curve edge drawings. We will confirm this 
equality after we describe the generation of curves in the next section. 

4.3 Curves 

After the positions of tree centers (and the positions of nodes if not given) are 
computed, we are now ready to place the confluent tracks for the edges. 

We use Bezier curves to draw the curve edges in confluent drawings. Given 
a set of control points Pi, P2, . . . , Pn, the corresponding Bezier curve is given by 

n 

=^PfcBfc,„(u) 0<M<1, (1) 
where Bk^n{u) is a Bernstein polynomial 

B^Au)^^^,u\l-uT-K (2) 

Bezier curves have some nice properties that are suitable for our confluent 
tracks. The first property is that a Bezier curve always passes its first and last 
control point. The second is that a Bezier curve always stays within the convex 
hull formed by its control points. In addition, the tangents of a Bezier curve at 
the endpoints are Pi — Pq and P„ — Pn-i - Thus it is easy to connect two Bezier 
curves while still maintaining the first order continuity: just let P„ = Pg 
the control points P„_i,P„ = PQ,and P{ co- linear. 

The confiuent track between each node and the tree center is drawn as a 
Bezier curve. In our program we use cubic Bezier curves {n = 4 in Equation 12) • 
Each such a curve has four control points, chosen as shown in Fig.^ 




Fig. 4. Bezier curves 

More formally, assume we are given the following input for a biclique P^: 
HmUi, and yc are the y-coordinates of the upper, lower, and tree center levels. 
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respectively. Xi is the x-coordinate of the tree center for Bi. the X- 

coordinates of nodes in biclique Bi . Let Ay be a distance parameter that controls 
the shape of the curve edges. When node j is in the upper level, the four control 
points are Pq = (xy , ?/„), Pi = {xij,yu + Ay), P2 = {xi,yc - Ay), and P3 = 
(xi,yc). When node j is in the lower level, the four control points are Pq = 
{xi,yc), Pi = {xi,yc + Ay), P2 = {xij,yi - Ay), and P3 = {xij,yi). 

From Equation^ it is not hard to verify that in a confluent layered drawing, 
two Bezier curves cross each other if and only if the corresponding straight-line 
edges (dashed lines in Fig. of the bicliques cross each other, given that the 
control points are chosen as above. This should clear the doubt that appears at 
the end of Section ^21 

5 Multi-depth Confluent Layered Drawings 

So far we have introduced the method of confluent layered drawings: replacing 
subsets of edges in a two-layer graph by tree-like structures. This method can 
be extended to obtain drawings that display richer information. The extended 
drawings are called multi-depth confluent layered drawings. 

The idea is as follows: after the biclique cover for a two-layer graph G = 
{U,L,E) is computed, the tree center nodes are viewed as a middle layer M, 
and a new three layer graph G' = {U, M, L, E') is constructed as in Section lOl 
The same biclique cover algorithm is then applied to G' twice, once for the 
subgraph induced by U \J M; once for the subgraph induced by M U L. By 
applying this approach recursively, we get biclique covers at different depth. 
In the final drawing, only biclique covers at the largest depth are replaced by 
sets of tree-like structures. The final drawing is a multi-depth confluent layered 
drawing. The drawings discussed before this section are all depth-one (confluent 
layered) drawings. 

In a depth-one drawing, we compute a biclique cover and lay out the biclique 
cover. In general, for a depth-i drawing, we need to compute 2^ — 1 biclique 
covers and biclique covers are laid out. 

An example drawing of depth-two is shown in Fig. [3 

Because the control points for the Bezier curves are chosen in a way such 
that the tangents at the endpoints of the Bezier curves are all vertical, it is 
guaranteed that all segments of a path are connected seamlessly and smoothly in 
multi-depth drawings. Readers probably have already noticed some wavy edges 
in the drawing of Fig.[3 It is because a single edge biclique is also drawn 

as two Bezier curves. We offer an option in our program to do a simple treatment 
for these single edge bicliques: draw them as a single Bezier curves instead of two. 
But after this special treatment is applied, the crossing property is not preserved 
any more. That means two curve segments could have crossing(s), even though 
their corresponding edges in G' don't cross each other in a straight-line drawing. 

Multi-depth drawings may further reduce the number of crossings. They also 
show a richer structure than the depth-one drawings, which only display bi- 
cliques. For example we can observe relationships between bicliques in a depth- 
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Fig. 5. Depth-two confluent drawing (on same input as the drawing of Fig. 

two confluent layered drawing. However higher depth requires more computa- 
tions of biclique covers, and generates more dummy centers. The former leads 
to the increasing of time and space complexity, while the latter could result in 
a more complicate confluent drawing. We feel that drawings with depth higher 
than two are not very practically useful. 

6 Real-world Examples 

We list two example drawings of real- world graphs in Fig.El We implemented the 
algorithm of computing biclique cover using the RLE heuristic. For the center 
placement we implemented the barycenter method. We assume that besides the 
two-layer graph, the input also includes the positions of (fixed) nodes in upper 
and lower levels (possibly output by other algorithms that take labels and other 
information into account.) The result drawing is written into a file of DOT 
format jl7| . The neato program in the Graphviz package [J] is then used to 
generate the graphic file in a desired format. Fig.|Hl(a) is a depth-one drawing. 
Fig. (b) is a depth-two drawing with the special smoothing treatment applied. 

7 Conclusions and Acknowledgments 

In this paper we introduce a new method - confluent layered drawing, for visual- 
izing layered graphs. It combines the layered drawing technique with the relaxed 
confluent drawing approach. There are still interesting open problems, e.g., how 
to test whether a layered graph has a crossing-free confluent layered drawing? 
How to minimize the crossing of the drawing among all possible biclique covers? 
It is also useful to investigate better ways for visualizing confluent tracks. 
We would like to thank anonymous referees for their helpful comments. 
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ciiids.h I /usr;'mcli]de/strmgs.h /usi/incliidc/ctypch /nsr/lncliide/stdiQ.h ctQQls.h 1 1 tcxchk.h tcxclLTis.h 




cmds.c ctQQls.c texclik.c verbatini.c cirors.c texchafs.c 

(b) "Includes" relation for the Texchk program 



Fig. 6. Confluent drawing for examples of Newbery |21j . 
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